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^ ^ Abstract. In this paper we consider an inverse problem for the time dependent linear Boltzmann equation. 

It concerns the identification of the coefficients via a finite number of measurements on the boundary. We 
CO , prove that the total extinction coefficient and the collision kernel can be uniquely determined by at most 

k measurements on the boundary, provided that these coefficients belong to a finite /c-dimensional vector 
space. 
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1. Introduction 



^ . In this paper we consider an inverse problem for the Unear Boltzmann equation 

(N 

OO ■ dtU + uj-W.^u + q{u-K^[u]) = in (0,T)xSxf7, (1.1) 



where T > 0, O is a smooth bounded convex domain of E^, A^ > 2, § denotes the unit sphere of M.^ , 
q € L^{i^) and K^, is the integral operator with kernel k{x,u!',uj) defined by 

KK,[u]{t,Ld,x) — / k{x,uj' ,uj)u{t,uj' ,x) duj' . 



}_] ' In applications, the equation (1.1) describes the dynamics of a monokinctic flow of particles in a body 

5^ I f2 under the assumption that the interaction between them is negligible. For instance, in the case of a 

low-density flux of neutrons (see [8], [11]), g > is the total extinction coefficient (the inverse mean free 

path) and the collision kernel k is given by 

k{x,ui' ,uj) = c{x)h{x,uj' ■ uj), 

where c corresponds to the within-group scattering probability and h describes the anisotropy of the 
scattering process. In this model, q{x)u{t,u!,x) describes the loss of particles at x in the direction lj at 
time t due to absorption or scattering and q{x)KK[u\{t,u!,x) represents the production of particles at x 
in the direction u from those coming from directions uj' . 

Our focus here is the inverse problem of recovery the coefficients in (1.1) via boundary measurements. 
More precisely, we are interested to recover q and k by giving the incoming flux of particles on the 



boundary and measuring the outgoing one. Since these operations are described mathematicany by 
the albedo operator Aq^K, (which wiU be defined in the sequel), a first general mathematical question 
concerning this inverse problem is to know if the knowledge of Aq^^ uniquely determines q, n, i.e., if the 
map (q, n) ^^ Aq^K. is invertible. 

We may precise this question. A first one is to know if the knowledge for Aq^ni^f) for all / determines 
(g, k) {inGnitcly many measurcracntcs); a second one is to know if the knowledge of Aq^nifj), for j = 
1, 2, . . . , fc, determines (g, k) (finite number of measurements) . 

There is a wide bibliography devoted to the first problem. We mention the general results obtained by 
Choulli and Stefanov [4]: they show that q and k arc uniquely determined by the albedo operator (see 
also [10]). We also mention the stability results obtained by Cipolatti, Motta and Robcrty (see [7] and 
the references therein). 

There is also a lot of papers concerning the stationary case (see for instance those by V.G. Romanov 
[12], [13], P. Stefanov and G. Uhlmann [14], Tamasan [15], J.N. Wang [16], and also the references therein). 

In this work we focus on the second question, concerning the identification by a finite number of 
measurements. Under certain hypothesis and assuming that K{t,ui\uj) = c{x)h{u}' ,u}), we prove that 
q and c can be uniquely determined by at most k measurements, provided that q and c belongs to a 
finite /c-dimensional vector subspace of C{Vt). To be more precise, we consider the initial-boundary value 
problem 

dtu(t,uj,x) + Lu ■ 'Vu{t,LU,x) + q{x)u{t,uj,x) ~ q{x)Ki^[u]{t,LU,x), 
<u(0,w,x) = 0, (cj,a;)G§xO, (1.2) 

^u{t,uj,a) :=^ f{t,uj,a), {uj,<T)eJ:-, t e (Q,T), 

where / is the incoming fiux and E^:= {{Lu,a) G § x dfl; w • J^(cr) < 0} is the incoming part of the 
boundary. Then our main result can be stated as follows: 

Theorem 1.1: Let ft C K^ be a bounded convex domain of class C^, M > 0, T > diam(O) and 
X: =■ spanjpi, p2, • • • : Pfc}; wfiere {/Oi, p2, ■ ■ ■ , Pk\ is a linearly independent subset of C(f7). We assume 
tliatqeL°°{il), \\q\\oo <M and k S L°°(f2; C(§ x S)). 

a) If q E X , then there exist cui, . . . ,u)k G § and fi, . . . , fk G Co((0, T) x S^j that determine q uniquely. 

b) If k{x,ui',u!) — c{x)h{uj' ,uj), where c G A" and h G C{§ x §) satisfies h{uj,Ljj) ^ for every u! G §, 
then there exist uji, . . . , 2^ G § and /i, . . . , fk G Co((0, T) x S^) that determine c uniquely. 

Remark: As we can see from the proof of Theorem 1.1, the functions fj,j~l,...,k, have the form 

„ . x _ j (f>j{<J- ia;)e'^(*~"-'^) in the case (a), 

•'■''^ ''^''^^■^ \(5<i^.(a)0j(a-tw)e'^(*-"-") in the case (b), 

where A > 0, t G (0,r), {uj,a) G S^, (f)j G C^(R^ \ il) and 5q. is the spherical atomic measure 
concentrated on ujj. The coefficients are identified by measuring the corresponding solutions on the 
outgoing part of the boundary, only in the directions wi, . . . , Wfc. 

The proof of Theorem 1.1 is based on the construction of highly oscillatory solutions introduced in [7] 
and some arguments already used by the author in [6]. 

We organize the paper as follows: in Section 2 we recall the standard functional framework in which 
the problem (1.2) is well posed and the albedo operator is defined; in Section 3, we introduce the highly 
oscillatory functions that will be used, in Section 4, to prove Theorem 1.1. 



2. The Functional Framework 

In this section wc introduce the notation and we recaU some weh known resuhs on the Transport Operator 
and the semigroup it generates (see [7] and the references therein for the proofs). 

Let ft C M.^ {N > 2) be a convex and bounded domain of class C^ and S the unit sphere of M^. We 
denote hy Q-. — Exft and E its boundary, i.e., E: = S x 90. For p G [1, +oo) we consider the space L'p{Q) 
with the usual norm 



l'"IU''(Q)-= ( / \u{uj,x)\^ dxduj\ 



where duj denotes the surface measure on § associated to the Lebesgue measure in R^ ^ . 
For each u G L'p{Q) we define Aqu by 

N Q 

{Aqu){uj,x):^ u ■ \7xu{uj,x) = ^<^k-^ — {^,x), lo = (cji, . . . , wat) 

fc=i ■^'^ 

where the derivatives are taken in the sense of distributions in O. 

One checks easily that setting 'Wp:= {u G LP{Q) ; Aqu G LP{Q)}, the operator (^Ao,Wp) is a closed 
densely defined operator and Wp with the graph norm is a Banach space. 

For every a G dfl, we denote t^(cr) the unit outward normal at cr G 90 and we consider the sets 
(respectively, the incoming and outgoing boundaries) 

E=^: = {(w, cr) G § X aO ; ± lj ■ i/(cr) > 0}. 

In order to well define the albedo operator as a trace operator on the outgoing boundary, we consider 
LP(E ; dS,), where d^: = \uj ■ i'{a)\daduj, and we introduce the spaces 

W^:={ueWp;u^^^ GiP(S±;dC)}, 

which arc Banach spaces if equipped with the norms 

i/p 



l"llw±- 



(\\urw,+I^J^-i^{T)Mu;,a)\Pdadco^ 

The next two lemmas concerne the continuity and surjectivity of the trace operators (see [2], [3] and [7]): 

7±:W±-^LP(ST;dO, l±{u):^ u\^^. (2.1) 

Lemma 2.1: Let 1 < p < +oo. Then there exists C > (depending only on p) such that 

\oj-iy{cr)\\u{uj,a)\Pdcrdoj<C\\u\f~^, Vu G W^^. (2.2) 

Moreover, if p > 1 and l/p+ 1/p' = 1, we have the Gauss identity 

divx(uvuj) dxduj — / oj ■ h'{(T)u{ij,a)v{uj,(T) daduj, (2-3) 

Q Jt, 



for all u e W^ and v G Wp. 

As an immediate consequence of Lemma 2.1, we can introduce the space 

)%:= {/e Wp; / \u; ■ i^{a)\\f{uj,a)f dudcT < +00} 
an we have that Wj^ = W^ = Wp with equivalent norms. 



Lemma 2.2: The trace operators 7± are surjective from Wp onto L^ {Y,^ ; d$,) . More precisely, for each 

p 

\\h\\^± < C'll/llLp(ST,dC), 



/ e LP{J^^;d^), there exists h e Wp such that j±{h) = / and 



where C > is independent of f. 

We consider the operator A : D{A) — > U'{Q), defined by {Au){lo, x): — ijj- \/u{lli, x), with D{A): = {u G 

Wp;«ls-=o}. 

Theorem 2.3: The operator A is m-accretive in LP{Q), for p G [1, +00). 

Corollary 2.4: Let f G L'p{Q), p G [1, +00) and assume that u G D{A) is a solution of u + Au — f. If 
/ > a.e. in Q, then w > a.e. in Q. In particular, it follows that 

I|m||li(q) < II/IIli(q)- 

It foUows from Theorem 2.3 and Corohary 2.4 that the operator A generates a positive semigroup 
{Uf){t)}t>o of contractions acting on LP{Q). 

Let q G L°°{il) and k:J7xSx§— >Mbea real measurable function satisfying 



\k{x,u}',uj)\ du}' < Ml a.e. Q x S, 
\k{x,uj' ,u!)\ du < M2 a.e. 51 x §. 



(2.4) 



Is 



Associated to these functions, we define the following continuous operators: 

1) B G C{LP{Q),LP{Q)) defined by B[u]{uj,x):^ q{x)u{u},x), 

2) Kk[u\{uj, x): — L k{x,uj' ,lu)u{lo' , x) duj' . 

It follows from (2.4) that K^ G C{LP (Q) , LP (Q)) Vp G [l,+oo) and (see [8]) 

WKMWlhq) < mI^"' m'^^^MlhQ) < max{Mi,M2}h|Up(Q). (2.5) 

The operator A + B — K^ : D{A) -^ LP{Q) generates a co-semigroup {U{t)}t>a on LP{Q) satisfying 

||(7(t)|U<e«, C:=||g-||oo+M2. 

We consider the initial-boundary value problem for the linear Boltzmann equation 

dtu(t,UJ,x) + Lo ■ 'Vu{t,LU,x) + q{x)u{t,uj,x) — qKK\u]{t,uj,x) 

m(0, u), x) = 0, (w, x) G S X O (2.6) 

u{t,uj,a)^ f{t,uj,(j), (cj,cr)GS", iG(0,r), 



where q G L'^{Q) and 

KK\u\{t,uj,x) ~ I k{x,uj' ,uj)u{t,uj',x) doj' , 



with K, satisfying (2.4). 

By the previous resuhs, it follows that, for / e Lp(0,T; LP{Yi^ , d,£_)), p E [1, +00), there exists a unique 
solution u e C{[0,T];Wp) H C^{[0,T];LP{Q)) of (2.6). This solution u allows us to define the albedo 
operator 

Ag^.:LP{0,T;LP{j:-,dO)^LP{0,T;LP{j:+,dO) 

Aq,K[f]{t,i^,cr):^u{t,uj,(T), {uj,<t) G E+. 

As a consequence of Lemmas 2.1 and 2.2, Aq^K is a linear and bounded operator. 

We also consider the following backward-boundary value problem, called the adjoint problem of (2.6): 

dtv{t,uj,x) + uj ■ Vv{t,u!,x) — q{x)v{t,uj, x) = —qK*[v]{t,uj,x) 

v{T, UJ, x) ==0, (w, x) G § X O (2.7) 

w(i,a;,cr) =5(i,w,cr), (a;,(T) G T.+ ,t G (0,r), 

where 5 G i^' (O, T; L^' (S+ , rf^) : P' ^ [1 , +00) and, 

K*^[v\{t,oj' ,x):— / k{x,uj' ,Lj)v{t,uj,x) duj, 
Js 

with the corresponding albedo operator A* ^ 

Al,:LP'{0,T;LP'{j:+,dO)^LP'{0,T;LP'{j:-,dO) 

•^q^K[9]{t,^,<^)-=^ v(t,uj,a), {uj,<t) G S~. 

The operators Aq_K and A* ^ satisfy the following property: 

Lemma 2.5: Let f G Lp{0,T; iP(S-;d^)) and g G Lp'{0,T; LP'(E+;d^)), where p,p' G (l,+oo) are 
such that l/p+ 1/p' = 1. Then, we have 



n(cj • v{(7))f{t, Lo, (T)A*q ^ [g] {t, UJ , (t) dadujdt 
_ J- 



n{uj ■ v{a))g(t,u!,a)Aq.n[f]{t,uj,(T)d(Tdujdt. 
J+ 

Proof: It is a direct consequence of Lemma 2.1. Let u{t,uj,x) the solution of (2.6) with boundary 
condition / and v{t, uj, x) the solution of (2.7) with boundary g. We obtain the result by using (2.3), once 
the equation in (2.6) is multiplied by v and integrated over Q. D 

As a direct consequence of Lemma 2.5, we have: 

Lemma 2.6: het T > 0, qi,q2 E L°°(f7) and ki,K2 satisfying (2.4). Assume that ui is the solution of 
(2.6) with coefficients qi,Ki and satisfying the boundary condition f G ^^(O, T; LP{T,~ , dS,)) , p G (1, +00) 



and that Uj ^'s the solution of (2.7), with q2,K2 and boundary condition g E L^ (0,T;LP (E+,(i^)), 
1/p + 1/p' = 1. Then we have 

n(q2{x) — qi{x))ui{t, CO, x)u2{t, co, x) dxdujdt 

+ / / {qi{x)Ki.,^[ui]{t,uj,x) - q2{x)KK,2[ui]{t,uj,x))u2{t,uj,x)dxdujdt 
Jo Jq 

n{uj ■ v{a)) [Aq^^m [f] - A2.K2 [/]] (*> ^, ^)9{ii w, a) dadudt. 

3. Highly Oscillatory Solutions 

In this section we present some technical resuhs related to special solutions of (2.6) and (2.7) that will be 
useful in the proof of Theorem 1.1. They were published in [7] and [5], but for the reader's convenience 
we present here their proofs. We denote by q the zero extension of q in the exterior of il, i.e., q{x) = 
for all X ^Vt. 

Proposition 3.1: Let T > 0, qi,q2 £ i°°(fi), and k satisfying (2.4). Let ■01, "02 e C{S,C^{R^)) such 
that 

supp-0i(cj,-)nn= (supp02(w,-) + 7'w)nII = 0, Vcj e §. (3.1) 

Tiien, there exists Co > such that, for each A > 0, there exist Ri,\ G C([0,r];>V2) and R*2 x ^ 
C([0,T];>V2) satisfying 

II^1,a||c([0,T];L2(q)) < Co, 11-^2, aIIc([0,T];L2(Q)) < Co, (3.2) 

for which the functions ui, Uj defined by 



ii(t, w, x): - 7Ai(cj, X - tuj)c- lo 9i(— -") ds^ix(t-u.^) ^ j^^^^^^ ^^ ^^ 
i* (i,w,x): = V'2(w, a; -tw)e-^o«^ ("-"") %-'^(*-"-)+i?*,^(i,w,x) 



(3.3) 



are solutions of (2.6) with q ~ qi and (2.7) with q ^ q2 respectively. Moreover, if k G L°°(il; L^(S x §)), 
then we have 

lim ||i?l,A|jc([0,T];L2(Q)) = lim ||i?2.Allc([0,T]L2(Q)) = 0. (3.4) 

A— 7-+OC A— > + oo 

Proof: Let u be the function 

u(i,a;,x): = Vi(^,a;-ta;)c"/o'«i(^-"")''"e'^(*-"-^)+i?(i,w,x). (3.5) 

By direct calculations, we easily verify that 

dtu + uj-\7u + qiu - qiKi^[u] = dtR + io ■ Vi? + qiR - qiK^[R] - e'^*giZi,A, 

where 

Zi,A(t, uj,x):^ I n(x, cj', w)V'i(a;', x - tu:')c' /o 9i(^-^'^') ''"c-'^^'-^dw'. (3.6) 



By choosing Ri,x e C^{[0,T];L^{Q)) nC{[0,T]; D{A)) the solution of 

' dtR + oj-\7R + qiR = qiK,[R] + c'^*(7iZi,A, 

i?(0, uj, x) = 0, (w, x) eS xn, (3.7) 

R{t,uj,a) = 0, (a;,cr)eS-, 

we see that (3.1) iniphes that the function u defined by (3.5) satisfies (2.6) with boundary condition 

Multiplying both sides of the equation in (3.7) by the complex conjugate of R, integrating it over Q 
and taking its real part, we get, from Lemma 2.1, 



-— / \R{t)\'^dLodx + - / uj-iy{a)\R{t)\'^dujda+ / qi\R{t)\'^dujdx 
'^dt Jq 2 Js+ Jq 

5R f qiK^[R]{t)R(t)dujdx = ^ 
Jo 



'Q 
It follows from the Cauchy-Schwarz inequality and (2.5) that 



e / qiZi^x{t)R{t)dujdx 



K4Rit)]\\R{t)\dxdLu < Ci||i?(t)||i2(Q), 
where Ci: = max{Mi, M2}. Therefore, we obtain 

^ll^(i)lli^(Q) < ^2||<Zl||oo||i?(i)||i2(Q) + |kl||oc||^l.A(t)||i2(Q), 

where C2: = 3 + 2Ci. Since i?(0) = 0, we get, by integrating this last inequality on [0,i], 

l|i?(t)lli.(Q)<lki||cocll^^"-^^^£||^i,A(r)|li.(Q)dr^ Vte[0,r]. (3.8) 

The first inequality in (3.2) follows easily because \Zi^x{t,i^,x)\ < ||'0i||coc"''^"°°"^Mi and, as the same 
arguments hold for 1*2 and i?2 a i "^^ ^^^^ obtain the second inequality. 

We assume now k e L°° (fJ; L^(§ x S)) . For each x G M.^ , the map uj' M> exp(iAaj' • x) converges weakly 
to zero in i^(S) when A — > +00 and the integral operator with kernel k(x, •, •) is compact in L^(§). So, 
we obtain from (3.6), 

lim ||Zi.A(t,-,a;)||L2(s) =0 a.e. in [0,T]xfi. 

Moreover, \\Zi^\{t, •, a;)||i2(g') < C, where C > is a constant that does not depend on A. The Lebesgue's 
Dominated Convergence Theorem implies that 

lim ||Zi,a||l2([o.t]xQ) =0. (3.9) 

A — ^+00 

From (3.9) and (3.8) we obtain (3.4), and our proof is complete. D 

7 



Corollary 3.2: Under the hypothesis of Proposition 3.1, ifqi,q2 G C(f7) and k G L°°(f7;C(S x S)), we 
have, for every cj G §, 

lim ||i?i,A(-,w,-)I|c([o,T];L2(t2)) = lim \\R2.xi-,^r)\\c{[a,T]:L^n)) =0. 

A— f+oo A— >+oo 

Proof: By multiplying both sides of the equation in (3.7) by the complex conjugate of R{t,uj,x), inte- 
grating it over 17, taking its real part and applying the Holder inequality, we get 



+ Ikilloo i\\K^[R]{t,uj)\\l,^,,^ + ||Zi,A(t,c^)||i2(o) 



(3.10) 



Since 



\K^,[R]{t,uj,x)\< I \k{x,J ,uj)\\R{t,uj' ,x)\du:' 
Js 

< ( \K{x,uj',uj)\duj'] ( / \K{x,uj',Lu)\\R{t,uj',x)\^ duj'] 

< Ml^'Mli^ (^JjR{t,uj',x)f duj'^ , 

we obtain 

||i^«[i?](t,c.)||i.(^) <Mi|[«:|U||i?(t)||i.(Q). (3.11) 

From (3.8), (3.10) and (3.11) we have 

^||i?(t,c.)||i.(,,)<4||gi|U||i?(t,c.)||i.(^) 

+ C (^||2'i,A|Ii2((o, T)xQ) + \\ZiAi't^^)\\h{n)) ■ 

Now, integrating this last inequality on time, we get 

||i?(i,^)|li2(o) <C7ell«^ll~^ (t\\Z,^x\\h^^o,T)xQ) + I ||^i,A(r,c.)||i.(,,)dT 

< C^IMUT (T||Zi,A||i.((o,T)xQ) + ll^l,A(-,a;, ■)lii2((o,T)xO) 

From Proposition 3.1 we know that ||.Z^i.aI[l2((o^t)xQ) ^ as A — > +oo. On the other hand, as the 
map uj' I— > e'" '^ converges weakly to zero in L^(S), we have from (3.6), for almost x E Q,, 

lim Zi^x{t,Lo,x) = 0, Vw G S, Vt G [0,T] 

A— 7-00 

and the conclusion follows from the Lebesgue's Theorem. D 

Lemma 3.3: We assume that q G L°°{V,) and k satisfies (2.4). Let S^ the solution of 

dtS + OJ-VS-qS= -qK*[S] + qe-'^'Z, 

S{T, UJ, x) := 0, (lu, x) G § X O, (3.12) 

Sit,u},a) = 0, (a;,cr)GS+, 



where Z e H^ (O, T; L'^{Q)) such that Z{T) = 0. Then we have 

I|5'a||c([o,t]:L2(q)) < Co and lim ||5^||h-i(o,T;L2(q)) = 0, (3.13) 

A— >oo 

where Cq is a constant independent of X. 

Proof: Multiplying both sides of the equation in (3.12) by the complex conjugate of S*^, integrating it 
over Q and taking its real part, we get 

~\\Sm\\UQ) + \ J^i^ ■ H<y))\Sl{t,LU,a)\duda > - \\q\\oo\\Siml.^Q^ 

- |k||oo|li^:[S](t)|U.(Q)||5I(t)|U.(Q) - |k|Ioo|I^(t)||L^(Q)|!^I(t)||L^(Q) 

Since \\K*[S]{t)\\L2(Q^ < niax{Mi, M2}||5^(t)|U2(Q), we have 
d 



dt 



WsmWhiQ) > -C2\\sm\\UQ) lkllcoll^(t)lli.(Q), 



where C2: = (3 + 2 niax{Mi, M2})||(7||oo. Integrating this last inequality on [t, T] and taking into account 
that Sl{T) = 0, we obtain 

WSimh^Q^ < hWoo^f^^ j^ ||^(r)||i.(Q)dT< ||g||ooe^^^||Z|U2(o,T;L^(Q)) (3.14) 

and the inequality in (3.13) follows easily. 
We consider now 

w\{t^u)^x):= I Sx{T,uj,x)dT, h{t,uj,x):= / c^' '^Z{t,oj,x) dr. (3.15) 

Then, it is easy to check that w\ satisfies 

dfW + Lu ■ \7w — qw = —qK*[w] + qh, 

w{T, Lj, x) ^ 0, (w, x) e § X rj, (3.16) 

w{t,uj,a) — 0, (w,o') G S+, 

Multiplying both sides of the equation in (3.16) by the complex conjugate of wa, integrating it over Q, 
taking its real part and applying the Cauchy-Schwarz inequality, we get as before, 

\\wx(t)\\l,^Q^ < ||g||coe^^^||/i||i.(o,T;L^(Q)) < \\q\\ooT'e^'^\\Z\\l,^,^^.j^.^Q)y (3.17) 

As S'{ = -dtwx, it follows from (3.14) and (3.17) that the set {wx} is bounded in C^{[0,T];L'^(Q)) 
and, in particular, is relatively compact in C([0,T]; L^(Q)). 

On the other hand, by integrating by parts the second integral in (3.15), it is easy to check that there 
exists C > (depending only on T) such that 

C 

II^IIl2(o, T:L2(Q)) < 1^11^11^1(0, T;L2(Q)). (3.18) 

Hence, by (3.17), it follows that ||ica||c([o.t];L2(q)) ^ as A — ;> cxd. Since the partial derivative in t, 
dt : C{\0,T];L^{Q)) -^ H^^(0,T; L^{Q)), is a continuous operator, there exists a constant C3 > such 
that 

II'S'a||h-i(0, T;L2(Q)) = ||9tWA||ff-i(0, T;L2(Q)) < C2 || Wa ||c(0,T;L2(q)) 

and we have the conclusion. D 



4. Identification by a Finite Number of Boundary Measurements 

In this section we assume that {pi,p2, ■ ■ ■ , Pk} is a given linearly independent set of functions of C{Q) 
and we denote X: = span.{pi, p2, . . . , Pk}- For each w e § wc consider Pq [pi] the X-ray transform of pi in 
the direction 2, i.e., 

/oo 
Pi{x + tuj)dt 
-oo 

and, for cache > 0, fie: = {x e M^ \II; dist(a:;, f7) < e}. 

The following Lemma, which the proof is given in [6], will be essential for the proof of Theorem 1.1: 

Lemma 4.1: For all e > 0, there exist ljj G § and (fij £ C^{Q,^), j — 1, . . . ,k, such that the matrix 
A = (aij), with entries defined by 

a2 



^D,[ft](x)0j(x)dx, (4.1) 

is invertible. 

Remark 4.2: It follows from Lemma 4.1 and the equivalence of norms in finite dimensional vector spaces 
that there exists C > such that, for all p E X, 



k 

Pcj,[p]{x)(f)^j{x)dx 



E 



R« 



>C\\p\\ 



In order to prove Theorem 1.1, we define, for < r < 1, the function x^ ^ § x § ~^ K as Xj.{uj,uj):^ 
P{rui,Ljj), where P is the Poisson kernel for -Bi(O), i.e., 

P{x,y): 



From the well known properties of P (see [9] ) , we have 



Xr{^,^) duj — 1, Vr e (0, 1), Vw e §, 

(4.2) 



lim / Xr{^i'^)4'i'^) duj = ip{uj), 

where the above limit is taken in the topology of LP{§>), p E [1, +oo) and uniformly on § if V-" G C(§). We 
are now in position to prove Theorem 1.1. 

Proof of Theorem 1.1: Let e: = (T - diam(r2))/2. 

Step a: We assume that ki = K2 = k and qi,q2 G X. We define ipi{uj,x) = ip{x) and, for 2 G §, 
V'2(w,x) = Xr(w, w)V'(a;), where < r < 1 and (fi^ip E Cl^{fle), with HV'IIl^ < 1- Then V'l and 'ip2 satisfy 
the condition (3.1) and we may consider the solutions ui and 1*2 defined by (3.3), i.e., 

ui(i,c^,x):-(^(x-tc^)e"/o'«^^("-™)%'^(*---")+i?i,;,(i,c^,x), 
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where A > will be chosen a posteriori. We shall write 

$i,A(t,^,x):=^(a;-tc^)e"/o«^^("-™)%'^(*---) 

^2 At, c., x): = 4>{x - tc.)e/o ^-'^---')d- ^-^Kt---^) 
in such a way that 

Substituting tii and 1*2 in the identity given in Lemma 2.6, we have 

/(A,r)-J(A,r)-L(A,r), (4.3) 

where 

I{X,r):= / / (^q2{x) — qi{x))ui{t,uj,x)u2{t,UJ,x)dxdujdt, 
Jo Jq 

J(A, r):= / / (^q2{x) — qi{x))Ki^[ui]{t,uj,x)u2{t,uj,x)dxdujdt, 
Jo Jq 

L{\,r):= [ [ {lj ■ i^{a)){Ailfx] - A2lfx])92.x.rdadijdt. 
Jo Jt.+ 
In the above formulas, we are denoting Ai = Aq-, i = 1,2 and 

/A(i,..,a): = ^(a-t^)c-/o'^^^(^-™)%'^(*— ), (c.,a) eS^, 
52.A,r(i,c.,a): = x,-(2,c.)^(a-ta;)e/o «^(^-™)'^-e-'^(*— ), (c.,a) e S+. 

Since 51 is convex, for s > and (a;,(T) G S^, we have {a — slu) ^ 17 and qi{(7 — slu) = 0. Hence f\ does 
not depend on gi, i.e., 

/A(t,cj,cr) =^(cr-ia;)e'^(*-'^-"), {u},a)eT.-. (4.4) 

By denoting /9(a;) ~ §2(2;) ^ 91(2;) and by considering the special form of ui and U2, we may write 
J(A, r) and J(A, r) as / = /i + /2 + /a + /4 and J = Ji + J2 + J3 + J4, where 

Ii(X,r):— / / p(x)c Jo Xr{'^,^)'p{x — tuj)ip{x — tuj) dxdwdt, 



/2(A,r) 
/3(A,r) 
/4(A,r) 



p(a:;)<i>i,A(^,<^,a;)i?2.A rl^j'^i^) dxdurdt, 



pix)Xri^ , ^)'^2,\it, (^ , x)Ri,\{t, uj , x) dxdujdt, 



10 JQ 
i-T 



JQ 



p{x)Ri^\{t, uj, x)R2 A r(^j "^j 2;) dxdujdt. 



and 



Ji(A,r):= / / p(x) 

JO JQ 

J2(A,r):= / /p(x) 

JO JQ 

U\,r):= f f p{x) 

Jo JQ 

U\,r):= f f p{x) 

JO JQ 



k{x, a;',a;)$i,A(i,'^', x) duj' 



k{x,uj' ,Lu)^i_x{t,Lu' , x) dw' 



k{x,u)' ,u))Ri,\{t, Lo' , x) dw' 



k(x,u) ,uj)Ri,\(t,io ,x)duj 



Xr{'^T'-^)'^2.\it, Lo, x) dxdujdt, 
i?2 A r{t,'^, x) dxdujdt, 
Xr{uj,ijj)'^2.\{t,'^iX) dxdujdt, 
i?2 A r(t, UJ, x) dxdujdt. 
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Taking the limit as r — > 1 in the above expressions, we get from (4.2), Ii{X, r) — > Ii{X) and Ji(A, r) 
Ji{\), i =■ 1, ... ,4, where 

/i(A):= / / p(£)c Jo '"■'^ '"^■' '^iflx — tu})ip{x ~tuj)dxdt, 

Jo JQ 

^2(A):= / / p{x)^i,x{t,Lj,x)S2 x{'t7^7^)d-xdujdt, 
Jo Jq 

hW-^ / / p(x)'^2,\{t,Lo,x)Ri,x{t,io,x)dxdt, 



10 Jn 

r-T 



/4(A): 



p{x)Ri^\{t^Lj^ x)S2 \{t, Lo, x) dxdiodt, 



JQ 



Ji(A):- / / p{x) 
Jo Jn 

J2(A):- / f p{x) 
Jo Jq 

^3 (A):- / f p{x) 
Jo Jn 

MX):= f f p{x) 
Jo Jq 

and 5*2 )^ is the unique sohition of 



k{x, uj' ,Lu)^i\{t,Lu' , x)duj' 

s 

k(x,lj ,uj)Ri^x{t,uj ,x)dw 
k{x, lu' , u))Ri,\{t, (x>', x) dw' 



^2,A(i,i^,a;) dxdt, 
S2 a(^i'^7 2;) dxdijodt, 
^2.\{t, ^, x) dxdt, 
S'2J^(t, cj, x) dxdivdt, 



(dtS + LJ-VS- q2{S - k:[S]) = e-'^*(72^2,A, 

S{T,Lo,x)^0, {Lu,x)eSxn, 
S{T,uj,a)^0, (w,cr)eS+, 



where 

Moreover, from (4.4) and (4.2), it foUows that L{\,r) -^ L{X), where 



Z*.A(t, cj, x): = k{x, lu, u)ij{x - tu)Jo i2ix-sa,)ds^ix^., 



L(A): 



Jan 



(2 • i^{a))+ {Ai [fx] - A2 [/a]) {t, 2, (T)*2,A(t, ^, <y) dadt, 



where »4i[/A] denotes the zero extension of ,44/a] on dH,. Therefore, by taking the hmit as 

(4.3), we have 

4 4 

^J.(A)-^J,(A)=L(A). 



So, 



|/i(A)|<^|/.(A)|+^|J,(A)| + |L(A)| 



i=2 



In what follows, we prove that there exists C(A) > 0, C(A) — > as A — ^ 00, such that 



(4.5) 



(4.6) 
1^ in 



(4.7) 



^|/,(A)|+^|J,(A)|<C(A)||/ 



i=2 



i=l 



(4.8) 
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We begin with terms of odd indices. Remembering that HV'lloo < 1, it is easy to see that 

\hW\ < C||p||oo||i?l,A(-,'^,-)||c(0,T;L2(n)), 

\JiW\ < C\\p\\oo I I \K^[^i,x]{t,i:i,x)\dxdt 
Jo Jn 

k3(A)| < CIIpIU / / \K4Ri^x]it,^,x)\dxdt, 



(4.9) 



Jn 
where C = C{T, M). From CoroUary 3.2 it follows that 

lim ||i?i,A(-,<^,-)||c(o.T;L2(n)) = 0. (4.10) 

Since cj' h-> c'^^" converges weakly to zero in L^(§) as A — ^ oo, we have $i^a(^, ',2^) ^^ a.e. t and x. 
Since X^ is a compact operator in i^(§), we have Kf^[^i_\]{t, u),x) ^ Q a.e t and x. So, by the Lebesgue 
Theorem, it follows that 

lim f f \K,[^i^x]{t,^,x)\dxdt^O. (4.11) 

^^°°Jo Jn 

Since 

\K^[Ri,x]{t,^,x)\ < [\K{x,u;\u)Ri,x{t,uj',x)\duj' < \\k\\oo j \Ri.x{t,uj\x)\d^' , 
Js Js 

it follows that 

iJsWl < C||p||oo||K||oo||-Ria||c(0,T;L2(Q)) 

and we have from Lemma 3.1 

lim ||i?l,A||c(O.T;L^(Q))=0. (4.12) 

A— >-oo 

On the other hand, we have for the even indices: 

|/2(A)| < C||p||oo / / \f{x -tuj)S2x{t,uj,x)\dxdujdt, 
Jo Jq 

\hW\ < C||p||oo||-Ria||L2(0,T;L2(Q))||S'2_AllL2(0.T;L2(Q)), 

|J2(A)| < CIIpIIoo / / \K,[<i>i,x]it,Lj,x)S;x{t,i^,x)\dxdujdt, 
Jo Jq 

|J4(A)| < CIIpIU / / \K^[R,,xKt,oj,x)S^xit,oj,x)\dxdujdt, 
Jo Jq 

From Proposition 3.1 and Lemma 3.3, it follows that 

lim \\Ri,x\\l^o,t,lHq))\\SIx\\l^o,T;L^q))=0- (4.14) 

A— ^oo 

Moreover, since 5*^ is bounded in L^{0,T;L^{Q)), K^^i^x] -^ and K^Ri.x] ^ in L'^{0,T:L'^{Q)), 
and we have 

lim / / \K4^i,x]{t,io,x)S2x{t,^,x)\dxdujdt = 0, 

V (4.15) 



(4.13) 



^^°°./0 JQ 



lim j j |i4r„[i?i_A](^,'^,a^)5'2 ;)^(t,cj, a;)| dxcfwdt = 0. 
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On the other hand, since f G C^{Vt^), it follows from the choice of e that the map (t,cj,x) ^^ (p{x — tuj) 
belongs to H^{0,T; L'^{Q)). Therefore, by Lemma 3.3, 

\f{x - tuj)S2^^{t,UJ,x)\dxdujdt < \\(p\\H}-,iO,T;L^Q))\\S2,x\\H-^(0,T;L^(Q)) ^0 (4.16) 

JQ 

as A ^ oo and we conclude from (4.9)-(4.16) that there exists C(A) satisfying (4.8). So, it follows from 
(4.7) that 



|/i(A)|<C(A)||p|| 
Now we remark that 



Jdn 



{Ai[h]-A2[h]){t,Sj,<j] 



dtda. 



\h{X)\ 



/o Jo 



p{x)c Jo ""^ ^ '^{x — tCj)'4i{x — toj)dxdt 



1 _ e^ Jo ^(2'+'"^)'^'' 



(p{y)'ip{y)dy 



(4.17) 



By taking the supremum on ip £ L°°(R ), ||V-'IIoo < 1, we get 



sup 



I _ Q- jg p(v+sai)ds 



v{y)^{y)dv 



I _ ^- Jg p(v+sai)ds 



\'p{y)\dy- 



Since 



^ p(y+suj)ds 



> 



p{y + suS)ds 



-MT 



we have from (4.17) (with Co = e" 

(■T 



-MT\ 



Co 



R"J0 



p{y + suj)\<fi{y)\ dsdy 



<C(A)|iH|c 



Jdn 



{Ai[h]-A2[h]){t,i:j,<j] 



dtda. 



(4.18) 



Since (suppcp + so;) n O = for all |s| > T, we have 



/ / p{y + s(^)\v{y)\ 

^R" Jo 



dsdy 



R« 



p{x)\ip{x — suj)\ dxds 
p{y + suj)\ifi{y)\ dsdy 

Pcj[p]{y)\f{y)\dy 



(4.19) 



From (4.19) and (4.17), we get 



Co 



Pci^[p\{y)\<fi{y)\ dy 



<C(A)||H|c 



/o Jdn 
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{Ai[h]-A2[h]){t,^,a) 



dtda. 



(4.20) 



We are now in position to conclude the proof of step a. First of all, we consider in (4.20) w as the 
directions uji, . . . ,ujk and (p as the functions (pl, . . . ,(p1 given by Lemma 4.1, in such a way that (see (4.4)) 

/j,A(t,c^,^) = </','(^-Me'^^*""'"^- 
Hence, by Lemma 4.1 and Remark 4.2, we can write, for some C > 0, 

C\\qi - q2\\oo < C{X)\\q^ - q2\\oo 



E/ / {Ai[fj,x]-A2[fj,x])it,i:;,,a) 



j^^ Jo Jdn 



dtda. 



Therefore, \i Ai[fj^\\{t,Uj^a) = ^2[/j",a](^j Wj7ct) on S^ , for j = l,...,/c, we have 

C||<Zi-rZ2||oo<C(A)||(7i-g2!|oo 
and the conclusion follows easily if we choose A > large enough. 

Step b: We assume that qi = q2 = q and Ki(x,Lo' ,uj) = Ci{x)h{Lu' ,lu), where ci,C2 G X. For cD G §, we 
define V'llwjx) = Xs(w,a;)(/)(x) and ■i/'2(w,x) = Xr('^, w)0(a;), where < r, s < 1 and G C^(i7g). Then 
'01 and "02 satisfy the condition (3.1) and we may consider the solutions ui and ^2 defined by (3.3), i.e., 

u, {t, LO,x):^ x.(S, c.)0(x - tLo)c- lo «-™)'^V^(*— ) + R^.^^^ (t, u, x), 
u;(t,cj,a;):=Xr(w,w)0(x-icj)e/o«("-™)%-i^(*---")+i?*_^,.(i,w,x), 
where A > will be chosen a posteriori. We shall write 

$A(t, c^, x): = <P{x - tcj)e- X' 9(™)'i-e-iA..^ 
*A(t, c^, x): = 0(x - i^)e/o 9(-™)'i-e'^^-" 



in such a way that 

ui{t,uj,x) = c'^*Xs('^,w)$A(i,w,x) +i?i,A,s(i,w,a;), 

U2(t,CJ,x) =C"'^*Xr('^,'^)*A(i,W,a;) + i?2,A,r(*>^:2;). 



(4.21) 



Substituting ui and U2 in the identity given in Lemma 2.6, we have 

J(A,r,s)=i(A,r,s), (4.22) 

where 

J(A, r, s):= / / q{x)(^ci{x) — C2{x))Kh[ui]{t,ui,x)u2{t,uj,x)dxdujdt, 
Jo Jq 

L{\,r,s):^ I I {lo ■v{cF))[Ai[f\^s\- A2[f\,s\)g\,Td(Tdujdt. 
Jo Js+ 

In the above formulas, we are denoting Ai ~ A^, i — 1,2 and 



fx.s{t,UJ,(y)--^ Xs(w,w)$A(t,'^,0-), ('^jCr) G S 

9\,r{t,uj,a):^ Xr(S,a;)*A(t,w,CT), (a;,cr) G S^ 
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(4.23) 



In particular, it follows from the definition of the Albedo Operator and (4.21), 
^i[/a,.] - A2[.fx,s] = Ri,x,s " R2,x,s, on (0,T) X S+ 



(4.24) 



By denoting 77(2;) — q{x)[Ci{x) — C2(x)j and by considering the special form of ui and 1*2, we may write 
J(A, r, s) as J = Ji + J2 + J3 + J4, where 



Ji(A,r,s):= / / ri{x 
Jo JQ 



h(oj\LLi)xs{uJ,^')^\{t,uj\x)duj' 
X Xri'^i'-^)'^ \{t,i-^, x) dxdujdt, 



J2(A,r,s):= / / -qix) 
Jo Jq 

J3(A,r,s):= / / ?7(x) 
Jo Jq 

J4(A,r,s):= / / 77(2;) 
"'0 Jq 



h{uj' ,u})xs{i^,i^')^x{t,uj' ,x)du}' 
h{ijj' ,uj)Ri,x,s{t, <^', x)duj' 
h{ijj' ,uj)Ri,x.s{t, w', x)duj' 



R2 X r(^: "^j 2;) dxdujdt, 
Xr('^, '^)^a(^, 1^, a;) dxdiodt, 
R2 xj.{t,uj,x) dxdcodt. 



Taking the limit as r — > 1 in the above expressions, we get from 4.1 Ji(A, r, s) — > Ji(A, s), where 



Ji(A,s):= / / 77(x) 
Jo Jo 

J2(A,s):= / ///(x) 
Jo Jq 

J3{X,s):^ f f rj{x) 
Jo Jfi 

J4(A,s):= / / ri{x) 
Jo Jq 

and 5*2 \ is the unique solution of 



h{uj',uj)xs (w, uj')^x{t, uj', x)duj' 



h{uj' ,Lj)xs{^, ^')'^x{t, ^' 1 x)duj' 



h{u)' , uj)Ri^x.s{t,'-tj' , x)duj' 



h{ijj' , uj)Ri^x,s{t, uj', x)duj' 



"^x{t,uj,x) dxdt, 
S2 xii^ "^j 2;) dxdiodt, 

"^xitji^jx) dxdt, 

S2 \{t, to, x) dxdcodt 



dtS + LO-VS-qS = -qX:^ [S] + C-'^'qZlx, 
S{T,Lo,x) = 0, {Lo,x)eSxn, 
S(T,io,a) = 0, (w,(7)gS+, 



Moreover, from (4.23) and (4.2), it follows that L{X,r,s) -^ L{X,s), where 



L{\,s):^ / (S-t.(a))+(A[/A,s]-^2[/A,s])(i,2,a)*A(t,S,a)dadt 
Jo Jan 

= / / (w-zy(CT))+(i?i,A,s(i,W,CT)-i?2,A,s(i,W,Cr))*A(i,W,Cr)(fCT(it, 

Jo Jan 



(4.25) 



(4.26) 



where Ai[fx.i] denotes the zero extension of Ai[fx,i\ on 9f2. Therefore, by taking the limit as r — > 1 in 
(4.22), we have 

Ji(A, s) + J2(A, s) + JsiX, s) + Ja{\, s) ^ L{\, s). 
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Now, it is time to take the limit as s ^ 1 . For the first two terms of the right hand side of the above 
identity, we get (for i = 1,2) Ji(A, s) — ;> Ji(A), where 



-^i(A):= / / ri{x)h{uj,uj)^xit,Lo,x)'^x{t,uj,x)dxdt 
Jo Jn 

= h{uj,uj) / / ri{x)(f>{x — tuj) dxdt, 
Jo Jn 

•/2(A): = / / rj{x)h{u)^u))^\{t,u),x)S2 \{t,io^x) dxdujdt. 
Jo Jq 



(4.27) 



On the other hand, the dependence on s in the other terms is given by i?i,A,s and R2,x,s, which are the 
solution of {j = 1,2) 

dtR + uj-VR + qR = qK^^ [R] + e'^*gZj- A,s, 

R{0, w, x) = 0, (cj, x) eSxn, (4.28) 

i?(0,a;,cr) = 0, (cj,o-)eS+, 
where 

Zj,x,s{t,i^,x):= / Kj{x,uj' ,uj)xs{(^,^')^x{t,co' ,x)duj' . 



(4.29) 

It is an immediate consequence of (4.2) and the Lebesgue's Theorem that, as s — > 1, Zj \ s — > Zj \ in 
C([0,r];L2(O)), where 

Zj^x{t,u!,x):^ Kj{x,ui,uj)^\{t,uj,x). (4.30) 

Hence, 

lim_ i?,. A,s = %A in C{[0,T]:L\Q)), 



where Sj^\ is the solution of 



dtS + uj-VS + qS = qK^^ [S] + e'^*gZ,- A, 
S'(0,a;,a;) = 0, (w,.T)e§xf7, 

S'(t,(x),o-) = 0, ((x),o-)eS^, 



(4.31) 



and Zj_\{t,u),x):^ Cj{x)h{uj,uj)^\{t,uj,x). Therefore, Ji{X, s) — ^ Ji(A), (i = 3,4) and L(X, s) — > L(X), 

where 

'■^ r r /■ 1 

^a(^, i^, x) dxdt, 



h{ijj' ,u])Si,\{t,ijj' , x)diu' 
h{u!', uj)Si^x{t, uj', x)dujj' 



M\):= / ry(a;) 
Jo Jn 

J4(A):= / f v{x) 
Jo Jq 

L(A):= / / (S-i/(a)) + (^i,A(i,w,a)-52.A(i,w,(T))*A(t,':5,a)d(Tdi 
Jo Jen 



s 



S2 x{t,^,x) dxdiodt. 



(4.32) 



and we obtain 
where 



ki(A)|<|J2(A)| + |J3(A)| + |J4(A)| + |L(A)|, 



|./2(A)| < ||ry||oo||/i||coc^^ / / \^{x-tu)S-^),{t,uj,x)\dxdLodt, 

Jo Jq 

|./3(A)| < h||oo||<^||ooC^^||if,[5i,A]|U^(0,T;L^(Q)), 

\Ji{>)\ < \\v\\'>o\\Kh[Si^x]\\L2(0,T;L^Q))\\Slx\\L2(0,T;L^Q)), 

\L{\)\<\mooC^^ f I (<I5-i^(a))+|5i,A(t,C},a)-^2,A(t,23,a)|dadt. 
Jo Jan 



(4.33) 



(4.34) 
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Since (f> G C^(f7e), it follows from the choice of e that the function {t,uj,x) h^ (f){x — tuj) belongs to 
Hq{0,T;L^{Q)) (as a constant function on uj). Hence, we have 

|-^2(A)| < |b||ooe^'^||0||Hi(o,T;L2(Q))||S'2^AllH-i(O,T;L2(Q)). 

On the other hand, from the weak convergence to zero in L^ (O, T; L^{Q)) of 5i.a, it follows that 



lim \\Kh[Si^x]\\L^(a,T;L^(Q))=0. 



(4.35) 



Hence, we have from (4.33)-(4.35) and Lemma 3.3, 

\Ji{X)\ — \h{uj,uj)\ / / 'q{x)(p{x — tu))'^ dxdt 



Jn 



(4.36) 



Jo Jon 



where C(A) ^ as A — ;> +oo. 

Since (supp0 + suj) n O = for all \s\ > T, we have 



Jn 



rj{x)(f){x — tuj) dxdt 



K" Jo 



ri{y + suj)4>{y)'^ dsdy 



p{y + suj)(j){yf dsdy 



(4.37) 



PM{y)'f>{yf 



dy 



and we get 

\h{uj,Lj)\ 



P^MivWy? 



dy < C{\)\\r,\\ 



C2 / / (cj • i^(CT)) \Si^\{t,uj,a) - S2,\it,uj,a)\ dadt 
Jo Jan 

We are now in position to conclude the proof. First of all, we consider in the above inequality the direc- 
tions wi, . . . , Wfc and the functions (j)i, . . . ,(j)k given by Lemma 4.1, in such a way that, from Remark 4.2, 
we can write 



Co|!ci-C2|U<C7(A)||ci-C2|U + 



+ C2V' / / (w • z/((t)) \Si^\{t,ujj,a) - S2,\{t,ujj,a)\ dadt, 
^^^ Jo Jdn 

for some constant Co > 0. If we denote by 

Uij{t,uj,a) ^ Xs{^j,^)^\{t,uj,x) + Ri^x^s{t,uj,x), J = 1,2, j == 1,. . .,fc, 
it follows from (4.2) that, as s — > 1~, u^ -^ uf,, where 



(4.38) 



,# 



Su^^^ + S^^x, i = 1,2, j = l,...,fc 



and Scj is the spherical atomic measure concentrated on luj . It is clear that 



# # 



— Si^x — 5*2, > 



Therefore, if ufAt,u}j,a) ~ uljAt^ujj^cj) on E^i^., for j = 1, . . . , fc, it follows from (4.38) that 

Co||ci - C2II00 < C{\)\\ci - C2II00 
and the conclusion follows easily if we choose A > small enough. 



D 
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